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We discuss the modal properties of the r-modes of relativistic superfluid neutron stars, taking 
account of the entrainment effects between superfluids. In this paper, the neutron stars are assumed 
to be filled with neutron and proton superfluids and the strength of the entrainment effects between 
the superfluids are represented by a single parameter rj. We find that the basic properties of the 
r-modes in a relativistic superfluid star are very similar to those found for a Newtonian superfluid 
\ star. The r-modes of a relativistic superfluid star are split into two families, ordinary fluid-like 

. r-modes (r°-mode) and superfluid-like r-modes (r°-mode). The two superfluids counter-move for 

' the r'^-modes, while they co-move for the r°-modes. For the r°-modes, the quantity k = a/Q -)- m is 

I almost independent of the entrainment parameter rj, where m and a are the azimuthal wave number 

and the oscillation frequency observed by an inertial observer at spatial infinity, respectively. For the 
' r°-modes, on the other hand, k almost linearly increases with increasing rj. It is also found that the 

radiation driven instability due to the r'^-modes is much weaker than that of the r°-modes because 
the matter current associated with the axial parity perturbations almost completely vanishes. 

in 
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, The r-mode instability P, 01 belongs to the secular instability attributable to the so-called CFS mechanism H, ^ |^ 
^0 ' that drives unstable various modes of oscillation in a rotating star. The r-modes excited in neutron stars are expected 
— to work as a mechanism of decelerating the spin velocity of neutron stars by emitting gravitational waves, which will 
^ ' be detectable by LIGO and other detectors ^ recent review, see, e.g., fllj . 

^ In old and cold neutron stars having the interior temperature below the superfluid transition temperatures T ~ 

?H . lO^K, neutrons in the inner crust and neutrons and protons in the core are believed to be in superfluid states [T^ . It 
is therefore likely that most of observed neutron stars, which have cooled down very quickly via emission of neutrinos 
I> . ^3 1 have a core with superfluids. One of the important effects of rotating superfluids in cold neutron stars is mutual 
' friction, which provides a strong dissipation mechanism for the oscillation modes 0^ . The mutual friction is produced 
' by scattering of normal fluid particles off the vortices in the rotating superfluids. Lindblom and Mendell 53 calculated 
' the mutual friction damping associated with the r-modes in neutron stars with a superfluid core and concluded that 
the mutual friction in the core is ineffective to suppress the r-mode instability. Recently, Lee and Yoshida 131 have 
studied the r-modes in neutron stars with the superfluid core using a different numerical method, conflrming most of 
the results obtained by Lee and Yoshida 16] have also shown that the r-modes are split into ordinary fluid- like 
r-modes and superfluid-like r-modes, as first suggested by Andersson and Comer _17], and that the instability caused 
by the superfluid-like r-modes is extremely weak and easily damped by dissipation processes in the interior. Very 
recently, Yoshida and Lee |I8| have calculated inertial modes of superfluid neutron stars, and shown that the inertial 
modes are also split into two families, namely, ordinary fluid-like inertial modes and superfluid-like inertial modes. 
It is quite reasonable since the r-modes belong to a subclass of the inertial modes. Yoshida and Lee [3 found that 
all the inertial modes, except for the ordinary fluid-like r-modes, are strongly damped by the dissipation due to the 
mutual friction unless the entrainment effects between superfluids are extremely weak. 

Most of the numerical studies on the r-modes and the inertial modes in superfluid neutron stars have been done 
within the framework of Newtonian dynamics. Since neutron stars are general relativistic objects having the typical 
relativistic factor M{R)/R of order of 10"^ with M{R) and R being the gravitational mass and radius of the star, 
respectively, it is necessary to extend the Newtonian analyses to general relativistic ones. In this paper, employing 



'Electronic address: yoshidaQfisica. ist.utl.pt] 
t Electronic address: |lee@astr. tohoku. ac.jp[ 



2 



the two-constituent model develope d by Cater and his co-workers 0, |^ , and the formalism devised by 
Lockitch, Friedman, and Andersson |23L l24l| for inertial modes of relativistic normal fluid stars, we derive a general 
relativistic formulation for inertial mode oscillations in superfluid neutron stars. In Sec. II we present the basic 
equations employed in this paper for the dynamics of relativistic superfluids, and in Sec. Ill numerical results are 
given, and Sec. IV is for discussions and conclusions. In this paper, we employ geometric units, given by c = G = 1, 
where c and G are the speed of light and the gravitational constant, respectively, and sign conventions used in [251 . 



II. FORMULATION 



A. Two-constituent formalism for superfluid dynamics 



We assume neutron stars are composed of superfluid neutrons and a mixture of superfluid protons and normal fluid 
electrons. We assume perfect charge neutrality between the protons and the electrons because the plasma frequency 
of the mixture is much higher than the oscillation frequency considered in this study . The electrons therefore 
co-move with the protons. In other words, we do not need equations for describing the dynamics of the electrons. 
In the following, we loosely call the mixture of the protons and electrons "proton" . To describe the dynamics of 
the superfluids in a neutron star, we employ the relativistic two-constituent formalism, which has been developed by 
Carter and his co-workers |19l l20ll2lL 12^ . The fundamental quantity of Carter's superfluid formalism is the so-called 
master function. Although we have several choices for the master function, we take, following Comer, Langlois, and 
Lin '22*1, the total thermodynamical energy density —A as the master function. The master function A is assumed to 
depend on the three scalars = —nau", = —paP", and x'^ = —pcji'^, where n" and p" are the conserved number 
density current of the neutrons and the protons, respectively. 

Once an explicit functional form of the master function is given, the energy-momentum tensor Tg is given by 

T; = *<5^-fp"X/3 + "V/3, (1) 
where the scalar denotes the generalized pressure, defined by 

^--A-nV^-p^Xa, (2) 
and the one-forms and Xa mean the chemical potential covectors, defined by 

Ma = ■BriQ + Apa , Xa = + Cpa , (3) 

where 

(9A 5A 9A 

^~~d^' ^~ d^' ^~ d^- ^> 

The one-forms and Xa are conjugate momenta to and p", respectively. The quantities p, = (— /^c^")"'^/^ and 
X — (— XctX")^^^ are interpreted as the chemical potentials of the neutrons and the protons, respectively. From 
equation ((JJ, we can see that the thermodynamical quantity A determines the strength of the entrainment effects 
between the two supcrfiuids. In other words, if we have A — Q, there is no entrainment effect. The system of the 
dynamical equations for the two superfluids is composed of two continuity equations, given by 



and two Euler equations, given by 



Van" = , VaP" = . (5) 



«"V[,M/3] = , p"V[,X/3] = , (6) 



where means the covariant derivative associated with the metric tensor, and the square brackets denote anti- 
symmetrised averaging. Since and are the conserved currents, it is convenient to introduce two time-like unit 
vectors u°' and ti" defined by 

n" = nw" , = pv" , (7) 



where n and p denote the number density of the neutrons and the protons, respectively. 
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B. Equilibrium configurations 

Equilibrium states of a slowly rotating star are described by the stationary axisymmetric spacetime, given by the 
line element 

= g^pdx°'dx'^ = -e^'^'^'^'Ut^ + e^^^-'Ur"^ + r^{d6^ + sin^ Odip^) - 2ujr^ sin^ Qdtd^ . (8) 

Here, we have taken account of the rotational effects up to the first order of the stellar rotation angular frequency. In 
this study, we assume that the neutron and proton superfluids in equilibrium are in the same rotational motion with 
the uniform rotation angular frequency Q. The four- velocity of the two superfluids is therefore given by 

li" = 7 {e + flip"') = e-''M(i" + flip"') , (9) 

where and tp" are the time-like and rotational Killing vectors of the spacetime. In virtue of equation the 
conjugate momentum one-forms in equilibrium reduce to 

Ha = {Bn + Ap)ua = ^J-{n,p)ua, 

Xa = iAn + Cp)ua=x{n,p)ua, (10) 

and the energy-momentum tensor of an equilibrium star is given by 

T^" = -A u^'up + ^q^ = -A u"up + {A + fm + xp) Qp , (H) 

where denotes the projection tensor associated with the four-velocity w", defined by q'^j — b'^ \- vF^u^. 
Substituting equations (|l()|l into Euler equations we obtain the hydrostatic equations given by 



1 d/i dv 1 dx dv 

p, dr dr ' X dr dr ^ 



(12) 



which leads to 



^=C, (13) 
X 

where C is an integral constant. In this study, we take C = 1, assuming that the two superfluids are in chemical 
equilibrium in the equilibrium states. The metric coefficients are determined from the Einstein equations, which are 
written as 

dv o^ fM . \ dM 



4Trr* I , — — = -Anr^A , (14) 



dr V / dr 



d'^vj /4 dv d\\ dvj A f dv d\\ , , 

-:^h^-rhr + :7rh^-o, (15) 



dr^ \r dr dr J dr r \dr dr 
where 

M{r)^^{l-e-'^^), w = n-uj. (16) 

In this paper, we neglect a normal fluid envelope (e.g., ^2^), and we require that the two superfluids have the 
common outer surface whose radius R is defined by equations n{R) = rrin and x(i?) = rrip, where to„ and rrip stand 
for the rest mass of the neutron and the proton, respectively. The surface boundary conditions for the equilibrium 
structure are given as follows: 

1 / 2M{R)\ Rdvj 
..-In 1 ^j, -+3^ = ^^, at r = R, (17) 



where M{R) stands for the gravitational mass of the star. The boundary conditions at the stellar center are the 
regularity condition for all the physical quantities. 
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C. Perturbation equations for the two-fluid model 



To describe fluid perturbations in a star, it is convenient to introduce two kinds of changes in physical quantities, 
called the Eulerian and the Lagrangian changes Q . An Eulerian change SQ is the difference between the quantities 
Q in the perturbed and the unperturbed states at a spacetime point. The relation between the Lagrangian changes 
and the Eulerian changes are given by 

AnQ = SQ + L^„Q, ApQ = SQ + L^^Q, (18) 

where Lk stands for the Lie derivative along the vector A:", and we have introduced the two distinct Lagrangian 
displacement vectors and to describe the perturbed motion of the neutron and proton superfluids because the 
two superfluids can move independently. The Lagrangian displacement is regarded as a vector that connects fluid 
elements in the unperturbed state to the corresponding elements in the perturbed state. 

We suppose that perturbations of any physical quantity can be expressed in terms of the Lagrangian displacements 
and , and the Eulerian changes in the metric ha/3 = Sgap- According to similar consideration to that in ^, we 
can write the Lagrangian changes of the fluid velocities of the neutron and proton in terms of the Lagrangian change 
in the metric: 

A„w" = -^u'^u'^u^Ar^g^, , Apw" = -^u'^u'^u'' Apg^, , (19) 

where 

AngaP = haf3 + ^ a£.n[3 + V/3^„q , Apgafi = hafi + ^ a^p/l + ^ P^pa ■ (20) 

Here, we have used the relation — in the equilibrium. In this study, we consider no transfusion between the 
neutrons and the protons. The particle numbers of the neutrons and the protons are therefore conserved separately. 
Then, the conservation equations are given by 

^nn 1 „^ App 1 

= ^Angap , = - Apgap . (21) 

n 2 p 2 

The Eulerian perturbations of the velocities and the number densities can be expressed in terms of , and hap 

as 

5v^ = + ^u^u^'u^h^, = 6v^ + , (22) 

5n = -'^q'^^Kp-q1V^{n£,^,), 5p = -^q'^'' h^p - q'^V ^ip^ ■ (23) 

According to and assuming — in the equilibrium, wc can write the Eulerian perturbations of the conjugate 
momenta fia and Xa as 

dfla = A* (^^Ua + hapU^ + -UaU'^u" hf^,^ + pA{5Va — SUa) + {(-4 — A)Sp + {B — B)6n] Ua , 
^Xa = X i^Va + hapu'^ + -UaU'^u'^h^^ + nA{5Ua ~ 5Va) + { (C - C)5p + {A — A)5n] Ua , 

where 



(24) 



A = 


dB 


-2n^ — 
d-n? 


-2p^^- 
dp^ 


OA 


B = 


dn? 




^dA 
^ dx^ ' 




c = 


2p2 




2dA 





(25) 
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Note that the second terms on the right-hand side of equations (|24ll represent the non-dissipative drag force between 
the two superfluids, which is proportional to both the function A and the velocity difference Sva — Sua. From equation 
(jnj, the perturbed Euler equations are given by 

q^LuS^if3 - q^di3{u''dnp) + fj.Su^{df3Ua - daUp) = , 

q^LuSxp - q%iu''Sxp)+xSv''idf3Uc.-dc.up)^0, (26) 

where da means the partial derivative. 

It is convenient to introduce the vorticity equation when pulsations in a rotating star are considered. Making use 
of the conjugate momenta as dynamical variables, the vorticity equations can be obtained as 

^«rfMa/3 = , LydXal3 = , (27) 

where dfj^afj and dxafs are the exterior differentiation of the one-forms fip and xtii defined by 

d^iap = ^a^J.^3 - dp^la , dXa/3 = daXp " dpXa ■ (28) 

Then, the Lagrangian variation of these equations can be derived straightforwardly: 

AnLud^la/3 = LuAnd^la/S = , ApLydXaff = L^ApdXafS = , (29) 

where Andfiaf3 and Apdxaf3 can be written in terms of the Lagrangian changes in the conjugate momenta as 

A„d^laf3 = dAnflal3 = 9aA„^^ - dpAnfia , ApdXafl = dApXa/S = d^ApXp - dpApXa ■ (30) 

Note that since the equations u^dAnfia/i = and d^^dAnfipy] — {v'^ dApXaf3 — and d^^dApXp-y] — 0) are satisfied, 
dAnfiap (dApXap) has only two independent components. 

The metric perturbations h^p are determined by the linearized Einstein equations, given by 

SC^ = SttST;^ , (31) 

where 5Gp denotes the linearized Einstein tensor. Here, ST^ means the linearized energy-momentum tensor, and is 
given by 

<5T^" = S-^S-^ + Sn^fip + 6p"xp + n"6fip + p'^Sxp , (32) 

where 

5-^ = i(n"A*^ + p''x^)Kp - - K^Xo ■ (33) 



D. Equations of state: Expanded master function 

In this study, we make use of the same master function as that introduced in which is generally given by 

oo 

This expansion of the master function may be justified because — np = in equilibrium and the deviation from 
x'^ —np = for perturbed states is the same order of the perturbations. Although Comer and Joynt have recently 
discussed in greater detail the entrainment effects in general relativistic superfluids, we use the master function given 
above for simplicity. In terms of the master function H34|) , we can obtain the thermodynamical functions in pulsation 
equations as follows: 

A=-Xdn\p^), B=-PA-'-'-^. C^-^A-'-'-^, 

n n on p p op 

A = A+^, B = B+^, C^C + ^, (35) 
onOp On^ Op'^ 
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where we have used the relation = np in the equilibrium. Note that A becomes equal to A if all the expansion 
coefficients are separable in n and p in the sense that A; = fi{n) + gi{p) for appropriate functions fi{'n) and gi{p). 
Using the chemical potentials in the equilibrium given by 

9Ao 5Ao 



we have 



A-A='^ = '^, B-B='J^, C-C='^. (37) 
op on on op 

and, assuming all the expansion coefficients A^ are separable in n and p, we obtain 

A = A, B-B=^, C-C = ^. (38) 
an dp 

E. Oscillation equations for inertial modes 

If we assume the star in the equilibrium is stationary and axisymmetric, the time and azimuthal dependence of the 
perturbations can be given by exp{iat + im(p), where a is the oscillation frequency observed by an inertial observer 
at spatial infinity, and m is the azimuthal wave number. Because of the rotation effects, it is generally impossible to 
achieve the separation of variables for the perturbed quantities. We expand the perturbations in terms of the tensor 
spherical harmonics with different Vs for a given m. In this paper, we consider the oscillation modes associated with 
TO > 2 because we are interested in the modes that are unstable against gravitational radiation reactions. Thus, 
we can select the so-called Regge- Wheeler gauge in order to fix the gauge freedom for the metric perturbations [29|. 
Then, the metric perturbations are expanded as 

e^-'Hoiir) Hu{r) \ 

sym sym sin^ 9Ki (r) / 

if^oi' ir)^d^Y^7'{e, if) ihw{r) sin edeYi7'{0, ip) \ 

^hu'{r)^d^Y^T{9,^) z/zu.(r) sin ^9,1^^(0, ^) L-* ^ (39) 
sym / 

where Yf"{9,(p) is the spherical harmonic function, and I =^ \m\ + 2{k — 1) and I' = I + 1 for even modes, and 
I — \m\ + 2k ~ 1 and I' = I — 1 for odd modes, where fc = l,2,3,---. The even and odd modes are characterized by 
symmetry and antisymmetry of the eigenfunctions with respect to the equatorial plane. The Lagrangian displacements 
^na and can be expanded in terms of the vector spherical harmonics as 

inn^r^t = , innCpt = , (40) 

^ p2A „2A 

iKnC„r= J2 —Wni{r)Yr{e.v)e'^\ iKn^p,= ^ —Wpi{r)YriO,^)e'''\ (41) 

/>|m| />|m| 




it^n^ne = V (Vni{r)dgYr{e, ^) - lUni'-^d^YiT'ie, ^)] e"^* , 

U^rn\ V Sine J 

E {yvi^r)deYr{9M-iUpv^.d^yir((i,p>)\e'^\ (42) 



, , V smf 

l,V>\m\ ^ 



inn^n^ = J2 {Vni{r)d^Yr{e,v)+iU^i,sinedeY,7'{0,^))e'''\ 

l,l'>\rn\ 

inVt^p^ = Yl {Vpi{r)d^yr{e,^)+iUpi.smedeYiV'{e,^))e'''\ (43) 



l,l'>\'rn\ 
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where Kfl = a + mfl is the oscillation frequency in the corotating frame of the star. Note that the gauge freedom 
for the Lagrangian displacement is fixed so as to satisfy = Cpt = 0- The Eulerian perturbations of the number 
densities for the neutrons and protons are given by 

5n^ J2 Sniir)Yrie,ip)e'''\ Sp = Spi{r)Yr{e,'f)e'''' . (44) 

i>|m| ;>|m| 

When we are interested in inertial modes in a slowly rotating superfluid star in the lowest order in the rotation 
frequency Q, it is rather easy to generalize the formulation devised by Lockitch et al. 23] for rclativistic inertial modes 
in normal fluid neutron stars. According to |23j| . we assume that the perturbations satisfy the following ordering laws 
in the slow rotation limit of — > 0: 

i/l - 0(1) , Wn = 0(1) , Wp = 0(1) , Vn = 0(1) , Vp = 0(1) , 

Ho = o{n) , H2 = o{n) , K = o{n) , Sn = o{n) , Sp = o{n) , (45) 



ho ^ Oil), Up = 0{l), Un = 0{l), (46) 

and K — 0(1). Comer fsoj showed that solutions obeying these ordering laws are allowed in the perturbation equations 
for rotating stars, exploring the so-called zero- frequency subspace of eigensolutions to perturbation equations in a non- 
rotating rclativistic superfluid star. If we consider no rotational effects, the solutions subject to the ordering laws 
(|45|l are interpreted as infinitely degenerate g-modes of zero-frequency, while the solutions subject to the ordering 
laws (|46|l are consider to be a rclativistic and superfluid counterpart of the trivial toroidal mode in a non-rotating 
Newtonian normal fluid star [Blj . 

Substituting the perturbed quantities into linearized equations (|23() , H29() , and (|31() and assuming the ordering laws 
for the eigenfunctions (|45|1 and H46I) . we can obtain a system of infinitely coupled ordinary differential equations for 
inertial modes in a slowly rotating superfluid star. To write down the oscillation equations for the inertial modes, it 
is convenient to use vector notation for the eigenfunctions w„, Wp, v„, Vj,, u„, Up, h, and H, whose components are 
given by 

W„.k = Wnl , Wp,k = Wpl , Vn.k = Vnl , Wp.fe = Vpl , Un,k = U„l' , Up,fc = Upl' , hk = hoi' , Hk = Hu , (47) 

where I = |r7i| + 2(fc — 1), T = ^ + 1 for even modes and / = \m\ + 2fc — 1, Z' = / — 1 for odd modes, where k — 1, 2, 3, • • •. 
The particle number conservation equations H23|) for the neutron and proton are then written as 

n„ d(e~^'^w„) / r dn dX ^ du\ 

e^^r— — = - 1 + — - +r— + 2r—] 



ar \ n ar dr ar J 

9„ d(e'^'^'Wj)) ( r dp dX dv 
o^^r^ . P£ = - IH L j^r— + 2r— 



e'-r^- ^ = - l + --i:+r— + 2r— Wp + AoVp . (48) 

dr \ p ar dr ar J 

Independent components of the vorticity equations (|29|l lead to 

Liu„ - qMoVn + qKoWn + (Up - u„) + h = , 

LiUp - qMoVp + qKoWp + ^2(u„ - Up) + h = , (49) 

dr dr dr \ dr J 

dq \ ( _ dq 

qKiKi - i"-^Mi 1 u„ + m I g - r—A^ 1 v„ - rH 



'-^(vp - v„) + ^1 I e (v^f„ - v^fp) 4- e ''^(e Vp - e a 



dr dr dr \ dr ' 



qKiAi - '"^^^i j Up + TO - r^A^ ^ j Vp - rH 
r^(v„-Vp)+^2 ('e^\wp-w„) + e^''r|-(e-^''v„-e-2-Vp) ) =0, (50) 
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where 

d 



The linearized Einstein equations are reduced to 

rH = -K^^lQ^Te^\'^{n^l^Nn +PX^p) , (52) 

d f d _ dv ^ d\\ d ( 2\k , 2\ dv d\\ \ 2, 



r— (^T-^hj - (^1 + ^d^^" [^e^^Ai + 2 - + ''^ + ''-^ j 1 j h = 167re^^r^(n/^u„ + pxUp) (53) 

where 1 stands for the unit matrix. Here, the matrices Lqj ^Ij Ki, Aq, Ai, Mq, and Mi are defined in Appendix. 
The oscillation equations we use are given as a set of linear ordinary differential equations for the variables w„, Wp, 
u„, Up, h, and r^, which are obtained by eliminating the vectors v„, Vp, and H in equations H48|) . H5U|I . and H53|) 
using the algebraic equations H49|) and H52|l . 

The boundary conditions imposed at the center of the star are regularity condition that all the perturbation 
functions are regular and do not diverge at r = 0. This implies that the eigenfunctions we solve must vanish at r = 0. 
In order to determine the boundary conditions at the stellar surface, on the other hand, we follow the arguments 
similar to those given in p^ . In our specific equilibrium models, for which we have assumed an equation of state 
that is separable in n and p, the generahzed pressure of the fluids can be written as ^(n,p) = '^n{n) + ^p(p), 
where ^n(n) and ^'p(p) represent the pressures due to the neutron and proton superfluids, respectively. Since the 
two superfluids can flow independently, as appropriate surface boundary conditions we require that the Lagrangian 
changes in each of the pressures vanish at the free surface of the star, which leads to 



A„*„-^A„n = 0, Ap*p = ^A„p = 0. (54) 



The conditions A„n = and App = at the surface therefore result in the surface boundary conditions given by 
Wni = and T^^p; = for inertial mode solutions in the lowest order of fi. In the exterior of the star, the only 
non-trivial metric perturbations /iq;' are determined by equation H53I) . which has two independent solutions, one is 
regular at spatial infinity, the other singular. The regular solution can be given as a series expansion {2,1 124| : 



hm' ^Yhv^A-] , (55) 



s=0 



R 



I'+s 



where 

~ (/- + .-2)!(r + . + l)!(2r + l)! ( 2M{R) Y^ 

- .!(Z'-2)!a' + l)!(2Z' + . + l)! (^"^ ) '° ' ^ ^ 

and /i;/ is an arbitrary constant. If we require that the interior solution h^y must be continuous with the exterior 
solution H55|) at the stellar surface because the spacetime must be regular everywhere, we obtain the boundary 
conditions for the metric perturbations at the stellar surface given by 



lim 



hw {R + e) r-^hai, (R-e)- r^h^' (R + e) hoi' {R - e) 
ar dr 



0. (57) 



For numerical computation, oscillation equations of a finite dimension are obtained by disregarding the terms with 
I larger than Imax in the expansions of the perturbations, where Imax is determined so that the eigenfrequency and 
the eigenfunctions are well converged as Imax increases. We solve the oscillation equations of a finite dimension as an 
eigenvalue problem with the scaled oscillation frequency k using a Henyey type relaxation method (see, e.g., [Slf'). 

III. NUMERICAL RESULTS 
A. Equilibrium models 

For the equilibrium structure with — np = 0, we use a generalized polytropic equation of state, introduced by 
Comer et al. |22,j, which is given by 

Ao(n^,_p^) = -rn„n - (T„n^" - rupp - app'^" , (58) 
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where (t„, Up, and Pp are constants. Note that this master function is separable in n and p. The generahzed 
pressure and chemical potentials in the equilibrium can then be written as 

* = CT„(/3„-lK"+ap(/3p-l)/^ 

X = mp + CTp/Jp/''-^ (59) 
Because chemical equilibrium is assumed in the equilibrium state, we can analytically write p in terms of n as 

g \ l/(/3p-l) 



where we have assumed nip — rrin- We confirm from equations (|59|l and l|()U|) that two superfluids in an equilibrium 
state given by equation H58|l have the common outer surface if /3„ and f3p satisfy the conditions /?„ > 1 and /3p > 1. 
In order to introduce the entrainment effects in the perturbed states with x"^ — np ^ 0, we include, according to |27j| . 
the expansion coefficient Ai, given by 

M = -V —7 ^ 7, (61) 

rripp + f](mnn + rupp) 



where 77 is a parameter whose appropriate range is considered to be 0.04 < 77 < 0.2 [27|, |32, |33| ■ In this paper, we call 
77 the entrainment parameter. Making use of the formulas obtained in the preceding sections, we can explicitly write 
the relevant thermodynamical coefficients as follows: 

A = A = r, 7""^^ ^. (62) 

mpp + ?7(TO„n + mpp) 

In order to avoid unnecessary singular behaviors of p and n at the stellar surface, we consider only the case of 
Pn = Pp- Thus, all the models we use are non-stratified models in the sense that n/p = constant inside the star. The 
constant parameters cr„, Up, „ . a nd /3p in Aq for the equilibrium equation of state (see Table I) are the same as those 
used for model 1 calculated in |22| . In this paper, we consider two equilibrium models, whose physical parameters are 
tabulated in Table I. The model I is almost Newtonian whose relativistic factor M{R)/R is M{R)/R = 10~^, while 
the model II is a relativistic one with the relativistic factor M{R)/R = 0.15, the value of which is similar to those for 
the models used in [23|. Although the model I is not appropriate as a physical model of neutron stars, we use it to 
compare with completely Newtonian calculations [l^ [la | . 



B. r-mode oscillations 



We calculated the r-modes of the two equilibrium models for a wide range of the entrainment parameter. It is 
important to note that the V = m fundamental r-modes, whose eigenfunction Unm and/or Upm is dominating and 
has no nodes in radial direction, are the only r-modes we find in this study. This situation is quite similar to that for 
the r-modes in a barotropic ordinary fluid star, in which neither the overtone r-modes with V = m nor the r-modes 
with V 7^ 771, are found 24, 34, 35, 36]. It is generally found that the r-modes in a relativistic superfluid star are split 
into two families, which we call ordinary fluid r-modes (r°-modes) and superfluid r-modes (r'^-modes), according to 
the notation of Lee and Yoshida |0 . Because of doubling of the dynamical degrees of freedom for the system of two 
superfluids, the mode splittin g of this kind also appears in other oscillation modes essentially associated with fluid 
motions [1^ [l3 lla 123 123 IS^j 113 • In Figures 1 and 2, we plot the scaled frequencies k, = a/Vl + m of the r-modes 
associated with m = 2 and 3 as a function of the entrainment parameter 77. As shown by the figures, although n of 
the r*-modes increases almost linearly with increasing 77, k, of the r°-modes is almost independent of -q. Comparing 
the two frequency curves for the models I and II, it is found that both the r°-mode and r'^-mode frequencies are 
strongly dependent on the relativistic factor M{R)/R of the equilibrium models. This is because the frequency of the 
r-modes is roughly proportional to an average of the function nj within the star, which reflects the general relativistic 
effects such as the rotational frame dragging and the gravitational redshift. Similar dependence of n on the parameter 
M{R)/R for the r-modes has been found for relativistic ordinary fluid stars 23, 24, 41, 42, 43, 44J. 

From Table II, where k's for the r°- and r*-modes associated with m = 2 and 3 are tabulated for several values 
of 77, we can confirm that k's for the r°-modes for model I are nearly equal to 2m/l'{l' + 1), the value found for the 
r-modes in Newtonian ordinary fluid stars. It is important to note that k's for the r°- and r'^-modes at 77 = are not 
equal to each other, and the difference between the k's increases as the relativistic factor M{R)/R is increased. In a 
Newtonian star, however, k's for the r°- and r*-modes at 77 = are equal to each other in the lowest order of fl, as 
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shown in The reason for the difference between the Newtonian and relativistic r-modes at 77 = is that even 
at 77 = there still exist couplings between two superfluid motions through the Einstein equations in the relativistic 
case, but no such couplings exist in the Newtonian case because all the Newtonian gravitational perturbations vanish 
in the lowest order of fl. 

In order to illustrate how the eigcnfunctions of the r°- and r^-modes in a relativistic superfluid star behave, it is 
convenient to introduce a new set of variables defined by 

„r nWni+p Wpi 

W+l = ■ , W-l = Wnl - Wpl , 

n + p 

nVni+p Vpi , . 

V+l = ; — , V-l = Vnl - Vpl , (63) 

n + p 

U+l' = — — — ^ — — , U-l' = UnV — UpV . 

n+p 

In Figures 3 and 4, we display four dominant coefficients U+2, U+4, and ho2 for the r°-mode, and U-2, W-3, 

U-4, and ho2 for the r'^-mode for the case of m = 2 and 77 = 0.1 for model II. Here, the normalization condition 
has been given by Unm = 1 B,t r — R. The coefficients that are not displayed in these figures have very small and 
negligible amplitudes. Figures 3 and 4 represent that the basic properties of the eigcnfunctions are quite similar to 
those of the r°- and r'^-modes in a Newtonian star with superfluidity [l^ll8|. The neutrons and protons co-move for 
the r°-modes, while they counter-move for the r'*-modes. It is also found in Figure 4 that the metric perturbation 
ho2 almost completely vanishes for the r^-modes, which is consistent with the fact that U+m ~ 0. This means that 
gravitational radiations due to r'^-mode oscillations are negligible at least in the lowest order in fi. It is important to 
note that the other coefficients W±3 and U±4 are not necessarily negligible compared with U±2 because of the general 
relativistic effects I2II2II4I. 



IV. DISCUSSIONS AND CONCLUSIONS 



In this paper we have studied the modal properties of the r-modes of relativistic superfluid neutron stars, taking 
account of the entrainment effects between the neutron and proton superfluids. To describe the general relativistic 
dynamics of the superfluids, we employed the two-constituent formalism developed by Carter and his co-workers 
[13. I20I I21I |2^ . We derived the perturbation equations for the relativistic inertial modes in neutron stars filled 
with the superfluids by generalizing the formalism employed by Lockitch, Andersson, and Friedman [24l] for 
the relativistic inertial modes in normal fluid neutron stars. We found that the basic properties of the r-modes in a 
relativistic star with the two superfluids are very similar to those in a Newtonian superfluid star [T^llSj . We confirmed 
that the r-modes of relativistic superfluid stars are split into two families, ordinary fluid-like r-modes (r°-mode) and 
superfluid-like r-modes (r'^-mode). The two superfluids counter-move for the r*-modes, while they co-move for the 
r°-modes. The dimensionless frequency k for the r°-modes is almost independent of the entrainment parameter 77. 
For the r'^-modes, on the other hand, n almost linearly increases with 77. The gravitational radiation driven instability 
due to the r^'-modes is much weaker than that of the r°-modes because the matter current associated with the axial 
parity perturbations vanish almost completely for the former. 

A solid crust near the surface of the neutron star has a significant infiuence on the modal properties of the r-modes in 
the superfluid core, since the solid crust supports its own oscillation modes EH E3j ISl j a-nd resonance phenomena 
are expected between the r-modes in the superfluid core and the torsional sound waves in the solid crust even in the 
general relativistic context jS^ [s^l ■ Dissipation in the viscous boundary layer at the interface between the fluid core 
and the solid crust is another important issue for the r-modc instability, as shown by '5 1*1 . How the dissipation in the 
core-crust interface can be important for the r-modc instability will be a quite interesting problem when superfluid 
neutrons and protons in the core and superfluid neutrons in the inner crust are taken into account simultaneously 
(see, e.g., HI). 

Newtonian superfluid neutron stars can support infinite number of inertial modes, which are als o sp lit into two 
families, ordinary fiuid-like inertial modes (i°-mode) and superfiuid-like inertial modes (i^'-mode) |l6l Il8| . Non-linear 
couplings between the r-modes and inertial modes will be important to limit the amplitude growth of the r-mode 
instability, as recently shown by Arras et al. [s^l for Newtonian normal fluid stars. To investigate the relativistic 
inertial modes in superfluid neutron stars will be one of our future studies. 
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Appendix: Matrices used in pulsation equations 

The components of the matrices, Lq, Li, Kq, Ki, Aq, Ai, Mq, and Mi are given as follows: 
For even modes, 

(T \ - 1 (r ^ -1 



+ ' (/ + i)(; + 2)' 



jm jm 



jm jm 

(Ao)i,, = «a + 1) , (Ai),,, = (; + i)(; + 2). 



I Z + 3 



where / = |to| + 2i - 2 for i = 1, 2, 3, 
For odd modes, 



(T \ -1 (T \ -^ 



l{l + l) ' ' ''''' 
jm 

(^o)m = — h> (^o)i 



"^i (T^ \ — __!±1. 



(i^i)i,. = -h > (^i)m 



(Ao)i,i = /(/ + 1) , (Ai)i,i = l{l - 1) , 



where Z = |m| + 2i - 1 for i = 1, 2, 3, • • •. Here, J/" is the function of m and Z, defined by J/" = [(Z^ - m^)l{Al'^ - 1)]^'"^, 
and g = 2vj/ kQ.. 
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TABLE 1: Parameters describing stellar models 1 and 11. 





model I 


model II 




0.2 


0.2 


(Jp/mn 


0.5 


0.5 




2.0 


2.0 


Pp 


2.0 


2.0 


Uc (fm"'^) 


0.0025 


0.672 


Pc (fm--^) 


0.001 


0.269 


M/Mq 


0.009 


1.033 


R (km) 


13.39 


10.17 


M{R)/R 


0.001 


0.150 



TABLE II: Scaled frequencies k for the r°- and r^-modes associated with m = 2 and 3. 



V 


model 


r°(m = 2) 


r"(m = 2) 


r°{m = 3) 


r"(m = 3) 





1 


0.6663 


0.6660 


0.4996 


0.4995 




II 


0.6094 


0.5459 


0.4411 


0.4194 


0.04 


I 


0.6663 


0.7592 


0.4996 


0.5694 




II 


0.6094 


0.6178 


0.4411 


0.4754 


0.1 


1 


0.6663 


0.8990 


0.4996 


0.6743 




II 


0.6094 


0.7240 


0.4411 


0.5585 


0.2 


1 


0.6663 


1.132 


0.4996 


0.8491 




II 


0.6094 


0.8964 


0.4411 


0.6942 



0.8 



1 1 1 1 1 1 1 1 1 1 

r°(model I) 


1 1 1 1 1 1 1 
m=2 mode 




. r*(model I) 






r''(model II) 






- r>(modeI II) 














, , 1 , 1 , , , , 1 







0.05 0.1 0.15 0.2 



FIG. 1: k's for r°- and r "-modes associated with m = 2 in the two models I and II, given as a function of tj. 
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r''(model I) 


m=3 mode 






_ r^(model I) 








r''(model II) 






- 


_ r^(model II) 




- 










1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



0.05 0.1 0.15 0.2 



FIG. 2: Same as Figure 1 but for modes associated with m = 3. 




FIG. 3: Four dominant expansion coefficients (7+2, M^+s, i7+4, and /io2 for the r°-mode associated with m = 2 in model II for 
r] = 0.1, given as a function of r/R. Other coefficients have negligible amplitude. The corresponding n is given by /t = 0.6094. 
The amplitudes are normalized by Un2 = 1 aX r = R. 
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FIG. 4: Four dominant expansion coefficients U-2, W-3, U-4, and ho2 for the r'^-mode associated with m = 2 in model II for 
q — 0.1, given as a function of r/R. Other coefficients have negligible amplitude. The corresponding k is given by /t = 0.7240. 
The amplitudes are normalized by Un2 = 1 at r = R. 



